function[EEF] = calc_eef(fday, 